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String breaking in Lattice QCD * 
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The separation of a heavy quark and antiquark pair leads to the formation of a tube of flux, or string, which 
should break in the presence of light quark-antiquark pairs. This expected zero temperature phenomenon has 
proven elusive in simulations of lattice QCD. We present simulation results that show that the string does break 
in the confining phase at nonzero temperature. 



In the absence of light quarks the heavy quark- 
antiquark potential is known quite accurately 
from numerical simulations of lattice quantum 
chromodynamics Q. At large separation R, the 
potential rises linearly, as expected in a confin- 
ing theory. In the presence of light quarks it 
is expected that the string between the heavy 
quark-antiquark pair breaks at large distance. All 
the existing lattice data at zero temperature pL[| 
agree in that they do not show any indication 
of string breaking which would be signalled by a 
tendency of the potential to level off at large dis- 
tances. The distances covered so far extend up 
to R<2 fm while it has been proposed that the 
dissociation threshold would be reached at sepa- 
rations somewhere between 1.5 and 1.8 fm 

We have simulated QCD with two light flavours 
of staggered dynamical quarks on lattices of size 
16 3 x 4 (new work) and 12 3 x 6 (configurations 
from Ref. |5)) at fixed values for the quark mass 
of m q /T — 0.15 and 0.075 respectively. The cou- 
plings were chosen to cover temperatures T below 
the critical temperature T c in the range of ap- 
proximately 0.7T C < T < T c . The (temperature- 
dependent) heavy quark potential V(R,T) was 
extracted from Polyakov loop correlations 



(L(6)L\R)) = c exp{~V(\R\,T)/T} 



(1) 



where L(x) = |tr hj^ 1 Uo(x, t) denotes the 
Polyakov loop at spatial coordinates x. In the 
limit R — > oo the correlation function should ap- 



* Poster presented at the XVI International Symposium on 
Lattice Field Theory, Boulder, CO, July 13-18, 1998. 



proach the cluster value |(L(0))| 2 which vanishes 
if the potential is rising at large distances (con- 
finement) and which acquires a small but finite 
value if the string breaks. 

In Figures |l| and || our data for the potential 
are presented in lattice units at the values of ft 
analyzed. The critical couplings (3 C have been 
determined as 5.306 for N T = 4 and 5.415 for 
N T — 6 respectively. The Polyakov loop corre- 
lations have been computed not only for on-axis 
separations but also for a couple of off-axis dis- 
tance vectors R. Rotational invariance is reason- 
ably well recovered if one uses the lattice Coulomb 
behaviour to determine the quark-antiquark sep- 
aration, \R\ = 1/GlatCR). 

The data in Figures |l| and || quite clearly show 
a flattening of the potential at lattice distances 
of about 3 to 4 lattice spacings, depending on 
ft. This confirms earlier results Q obtained on 
smaller lattices of size 8 3 x 4. Moreover, the 
height of the potential at these distances is in 
nice agreement already with the infinite distance 
cluster value, shown as the right-most data point 
in each of the plots. 

In order to obtain a rough estimate of the 
corresponding temperatures in units of the crit- 
ical temperature we applied the following proce- 
dure: at the given (3 and m q a values an inter- 
polation formula was utilized to estimate the 
vector meson mass mya in lattice units as well as 
the ratio of pseudoscalar to vector meson mass, 
mps/mv- By means of a phenomenological for- 
mula which interpolates between the (experimen- 
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Figure 1. The potentials in lattice units at the (3 
values analyzed for N T = 4. The right-most data 
points plotted at R/a = 9.5 and denoted by stars 
are the infinite distance cluster values — Tln| (L}\ 2 . 



tally measured) p and K* mass as function of 
the ratio m PS /mv, m v = 756MeV + 450MeV x 
i m PS l m v) 2 , a physical value for my can be ob- 
tained. These numbers are then used to estimate 
the ratios of the lattice spacing at the various (3. 

Finally, in order to facilitate a comparison of 
the N T = 4 and 6 results with each other and 
with quenched data, the absolute scale was de- 
termined from a conventional Wilson loop mea- 
surement of the string tension at zero tempera- 
ture at the critical (3 C values. The Wilson loops 
did not show string breaking at the separations 
which could be explored. The results for the crit- 
ical temperature in units of the string tension are 
obtained as T c j^f8 = 0.436(8) for N T = 4 and 
T c /y/a = 0.462(9) for N T = 6 g. 

In Figure [| we show the potential in the pres- 
ence of dynamical quarks in physical units. The 
potential is flat within the error bars at distances 
larger than about 1 fm. It also seems that the 
turn-over point is slightly T dependent, becom- 
ing smaller with increasing temperature. 

In Figure || quenched and full QCD potentials 
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Figure 2. Same as Figure 1 except for N T = 6. 



are compared. The quenched data has been taken 
from ||] and was obtained in the same way, i.e. 
computed from Polyakov loop correlations. Fig- 
ure |] contains, for further comparison, the dashed 
line denoting -n/(l2R) + (420MeV) 2 i? which 
gives a good description of the zero temperature 
quenched potential. Note that the finite temper- 
ature quenched potential is rising with distance 
R but the slope decreases with temperature, i.e. 
the (quenched) string tension is temperature de- 
pendent and becomes smaller closer to the critical 
T c . Again, the comparison with quenched poten- 
tials at the same temperature demonstrates that 
the potential in the presence of dynamical quarks 
becomes flat within the error bars at distances 
of about 1 fm. From Figure [| we conclude that 
the observed string breaking, albeit at finite tem- 
perature, is an effect caused by the presence of 
dynamical fermions. 

We have seen that string breaking is relatively 
easy to observe in the Polyakov loop correlation, 
while it is difficult to detect through the con- 
ventional Wilson loop observable. Why is this 
so? The Wilson loop observable creates a static 
quark-antiquark pair together with a flux tube 
joining them. In the presence of such a static pair 
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Figure 3. The potential in physical units at var- 
ious temperatures. The results are from lattices 
with N T = 4 and 6, indicated by the number in 
brackets. The data has been normalized to the 
cluster value. 



at large R, we expect the correct ground state of 
the Hamiltonian to consist of two isolated heavy- 
light mesons, however. Such a state with an ex- 
tra light dynamical quark pair has poor overlap 
with the flux-tube state, so it is presumably re- 
vealed only after evolution to a very large tem- 
poral separation. An improved Wilson-loop-style 
determination of the heavy quark potential in full 
QCD would employ a variational superposition of 
the flux-tube and two- heavy-meson states ■ 
The Polyakov loop approach, on the other hand, 
although limited in practical application to tem- 
peratures close to or above T c , builds in no preju- 
dices about the structure of the static-pair ground 
state wave function. Screening from light quarks 
in the thermal ensemble occurs readily. 
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Figure 4. The potential in physical units at vari- 
ous temperatures. Compared are quenched (open 
symbols) and full (filled symbols) QCD potentials 
at the same temperature. The dashed line is the 
zero temperature quenched potential. The data 
has slightly been shifted as to agree at distances 
around 0.3 fm. 
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